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Abstract. 

The odderon intercept a is recalculated by the variational method based on conformal 
invariance. The final result is a = 1 — 0.339 a s N/ir in contrast to the published value above 
unity. 
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1 Introduction 

As shown in [1,2], at least in the high-colour limit, all states formed by any number of 
reggeized gluons reduce to two basic ones: the pomeron for the positive signature and the 
odderon for the negative signature. The intercepts of these two states are fundamental for the 
high-energy behavior of cross-sections in the perturbative QCD. Unlike the pomeron inter- 
cept, known since more than two decades ago, the odderon intercept has not been calculated 
exactly until now, despite much activity around this problem and most sophisticated mathe- 
matical tools used [3-6] . The only published results refer to estimates made in the variational 
approach [7-9]. In paper [8], based on the conformal property of the odderon [7], the obtained 
intercept lies above unity. In the following we present it in the form 

a = l-— — -e (1) 

7T I 

Here a s is the strong coupling constant, N is the number of colours; e has a meaning of the 
odderon "energy" per gluon in units a s N/ir. The corresponding energy for the pomeron is 
—4 In 2. It is negative, which implies the intercept above unity and cross-sections growing 
with energy like s A where A = —(a s N /n)e. For the odderon in [8] it was obtained 

e = -0.25 (2) 

which is also negative but considerably greater than the pomeron energy. If correct, it would 
correspond to the negative signature cross-sections also rising with energy, although much 
slowlier than the positive signature ones. 

In the subsequent paper [9] the conformal properties of the odderon were not taken into 
account. Instead its wave function was presented as a sum of one-gluon "orbitals", which 
made it possible to take a very large number of variational parameters. A positive value was 
obtained for the energy: 

e = 0.298 (3) 

In principle, due to the variational nature of the calculations, this value does not contradict 
(2). However in [9] certain theoretical arguments were also presented in favour of the odderon 
intercept being smaller than unity. 

Guided by the evident conflict of these theoretical arguments with the result (2), we 
undertook to repeat the variatuional calculations in the conformal approach of [7,8], taking 
both the same trial function and the ones with more parameters. It turned out that the 
result (2) is in error. Correct calculations with the same trial function (with 2 parameters) 
give 

e = 0.2286 ±0.0001 (4) 
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with an opposite sign as compared to (2). Inclusion of more parameters (up to 8) allows to 
slightly lower this value down to 

e = 0.2262 ± 0.0001 (5) 

The uncertainties in (4) and (5) refer to errors in numerical calculations. In view of a very 
small change in e upon introducing more parameters, one is inclined to take (5) as a final 
value for the odderon energy with a precison of 0.1%. 

The study of the odderon energy in the conformal approach, in fact, presents a formidable 
calculational problem, deserving some attention. In this note we briefly discuss the difficulties 
involved and some intermediate steps in obtaining our results. 

2 Basic equations 

Basic equations for the odderon energy in the conformal approach were derived in [7]. Wc 
reproduce them here to fix our notatons. The odderon energy is to be sought as a minimum 
of the energy functional 

°° poo 

E= £ / dve n {v)\a n (v)\ 2 (6) 

Here 

poo p2n 

a n {u) = / drr' 2 ' 2iv / #e~^ 
Jo Jo 

(iu + + re»(h -iu- ^)) (iu - ^){-h + iv- ^)Z(r, <f>) (7) 

Z is the odderon wave function and 

e n {u) = 2 Re (yK^^ + ^) " ^(1)) (8) 
The normalization condition is 

°° POO 

1 = D= J2 / dv\a n {v)\ 2 (9) 

n=-oo J -°° 

In these formulas h and h are the two conformal weights, which are taken to be equal to 
1/2. The trial function Z then takes the form 

Z(r, (f>) = r 1 / 3 r 1 1/3 z(r, </>) (10) 

where 

r 2 = 1 + r 2 - 2rcos<f> (11) 

and z(r,(f>) = z(r,r*) shoud be invariant under the transformations of complex r: r — > — 1/r 
and r — > 1 — r. This can be achieved by taking z = z(a) where 

(12) 



(l+r 2 )(l+r2)(r 2 + r 2 ) 
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We take the trial function z(a) in the form 

N-l 

z( a ) = c k a k/2 ~ 1/6 + c N a 5/6 In a (13) 
k=i 

with N— 1 variational parameters ( one of is determined by the normalization condition). 
For N = 3 this function coincides with the one used in [8]. 

The basic quantity a given by Eq. (7) can be presented in the form 

a n {u) = q - «0(^ - + WPV) + (1 - \ ~ WfiMl (14) 

where 

/« = / drr- 1 - k - 2iu d^e- in(p Z(r,<f>) (15) 
Jo Jo 

Function fit\v) has the following properties. Since Z is real, (fn\v))* = f-l(— v) and 
since Z only depends on through cos0, we have fn\v) = f^(u). Under r — > 1/r the 
argument a stays invariant, from which it follows that Z — > Z/r. As a result we obtain 
fn\v) = fn\—v) = {fn^)*- Using these properties we can restrict the summation over n 
and integration over v to nonnegative values. The appearing coefficient 2 cancels in dividing 
E by D so that we need not take it into account. The value of the |a n (^)| 2 can evidently be 
expressed via a single function fn\v), which will simply be denoted as f n (v) in the following. 
In this manner we obtain our basic formulas for E and D: 

°° POO 

E = J2 due n {u)p n {u) (16) 

where for n > 

vniy) = 4 + v 2 ){^^ + + , 2 )|/ n (,)| 2 + (^^ + , 2 )|/„-iHI 2 + 

71 — I— 1 2 — 77 

2Re (-1- - «/)(— " "O/nM/n-lM] (17) 

and 

^^(I + ^I + ^i/^ (18) 

The normalization functional -D has the same form (16) with e n (v) — > 1. Thus calculation of 
the odderon energy requires calculation of functions / n (^) and e„(f). 

The main technical difficulty is a double Fourier transform (15). The energy e n {u) in 
Eq. (16), monotonously grows both with n and v. It is negative only for n = and small 
enough values of v. So the problem with this formalism is that cutting in (16) summation 
over n and integration over v by some maximal values fi m and u m , one always gets smaller 
E than the exact value, corresponding to n m and v m — > oo. Therfore in the course of the 
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calculation one always approaches to the variational value of e from below. As we shall see in 
the following, in fact, rather high values of n m and v m are necessary to obtain e with some 
degree of accuracy []. On the other hand, with high n and v, the double Fourier transform 
(15) becomes very difficult, especially having in mind that, due to the factors in (17), two 
first terms in the asymptotic expansion of a at high n and v cancel. As a result, as mentioned 
in the Introduction, a trustworthy calculation of E and D turns out to be very complicated, 
in spite of it superficial transparency. 

The crucial point in obtaining reasonable results has been using analytic asymptotic 
expansions for / at high n and v, which are discussed in the next section. 

3 Asymptotics at large n and v 

Passing to variable p = — lnr and introducing 2-dimensional vectors x = (p,4>). and w = 
(z, n) = (2v, n) we rewrite (15) as 

/ + OO /"+7T 
dp / dcf>e iwx Z(x) (19) 
-OO J — TV 

Evidently for the asymptotics at high n and v the integration point x = is essential. 
At x — > 0, keeping terms up to third order in small p and <f>, we have 

r 1 ^ = i _ (i/3) p + (l/18)p 2 - (l/162)/> 3 

rl /3 = x l ^(l - (l/6)p + (1/36)^ - (l/72)<f - (1/324)^ + (l/432)^ 2 ) 

a p = (1/2)V P (1 -p(29/12)p 2 -p(25/12)</» 2 ) 

(no terms of the third order appear in a). In these formulas x = \J p 2 + (ft 2 . As a result, at 
small x 

Z p = {mf'^aP = (l/2)V p+1/3 (l - {l/2)p + ap 2 + bcft 2 + cp 3 + dpcft 2 ) (20) 

where 

a = 5/36-(29/12)p; b = -1/72 - (25/12)p; c = -1/36 + (29/24)p; d = 1/144 + (25/24)p 

For the term with a logarithm, in the same manner we obtain 

Z p = (rr^WaP \na = (l/2) p x 2p+1 / 3 [ln(x 2 )(l - (l/2)p + ap 2 + bcft 2 + cp 3 + dp4> 2 )+ 

ln(l/2) - (1/2) ln(l/2) p + a lP 2 + hcb 2 + c lP 3 + d lP (ft 2 )] (21) 
1 This is the reason why in [7] the negative value (2) was obtained. The authors chose too small values of 
n m and v m 
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where 

h =/ln(l/2)+<5 / , f = a,b,c,d 

and 

S a = -29/12, 5 b = -25/12, <5 C = 29/24, <5 d = 25/24 

Inserting these expressions into the integral (19) and extending the integration over 4> to 
the whole real axis one obtains the asymptotical expansion of different terms in f n (v). In 
particular the asymptotical expansion of the term originating from Z v is found as 

f(p)(u) = c p (l + (l/2)id/dz - ad 2 /dz 2 - bd 2 /dn 2 + icd 3 /dz 3 + idd 3 / 'dzdn 2 ) W - a - 1 

where a = 2p + 4/3 and 

Cp = 2 a+1 - p T 2 {l/2 + a/2) cos(vra/2) 

Doing the derivatives, one obtains finally 

Re/#V) = c p w~ a - 1 (l + (a + l)ur 4 [z 2 (6 - a(a + 2) + n 2 {a - b(a + 2))]) (22) 

Iinf^V) = -(l/2)c p (a+l)^- a - 3 (l + 2{a + 3)u;~ 4 [z 2 (c(a + 2) - d) + n 2 (d(a + 4) - 3c)]) 

(23) 

For the term with a logarithm at p = 5/6 only the part with In a; 2 contributes. The 
result coincides with the formula above with p = 5/6 and the constant c p substituted by 
16vr(l/2) 5 / 6 . 

If one puts these asymptotic expressions into (17) one finds that the two leading terms 
coming from Z 1 / s cancel. Numerically the asymptotic expansion begins to work at rather 
high values of n and v: \Jn 2 + Av 2 > 30. 

4 Numerical procedure 

According to (15) function f n (u) is defined as an integral 

/•(X) r2n 

fn= drr- 1 - 2 ™ dcj>e- in( l>Z{r,4>) (24) 
Jo Jo 

The double Fourier transform (24) was done numerically in the interval 0<n<30, < ^ < 
15. 

At r > 1 we transform the integration variable r — > 1/r, under which Z — > Z/r. Then 
the integration over r becomes restricted to the interval < r < 1. Using properties of Z we 
can also integrate only over < < 7r and substitute e~ m< ^ by 2 cos n<p. We finally use that 
Z linearly goes to zero as r — > and introduce a new function 

C(r,0) = Z(r,$/r (25) 
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In terms of C we then have 



i 



f n = 2 I dr{r- 2iv + r ~ l+2iv ) I #cos in<f> ((r, <j>) (26) 



o 



The integrand is evidently singular at r = 0. To soften this behaviour we use that at 
small r 

C(r,(p) = Co-Cirlnr + C 2 r + 0(r 2 lnr) (27) 



where 



Co = c 1 (l/2) 1 /3 ; Cl = - 2cN (l/2f/ 6 , ( 2 = (l/2) 5/6 (c 2 - c N ln2) (28) 



In (26) we subtract from ( the first term of expansion (27) in the part of the integrand with 
r~ 2 ' w and the first three terms of this expansion in the part with r~ 1+2tu , adding to f n (v) 
the result of the integration of subtracted terms, which can be found explicitly. In this way 
we come to our final formula 

/» = 2 r#cosm</> f 1 dr (r- 2w (((r, 4) -(a) 
Jo Jo 

+r- 1+2 ^(C(r, <f>) - Co + Cir lnr - C 2 r)) + 

Wf( MT^j + (i + C L)2 + irk (29) 

Eq. (29) was used for numerical calculation of fniy) in the above mentioned interval of n 
and v. Integrations were performed by dividing the rectangle < r < 1, < < 7r into an 
M x M grid, interpolating C quadratically on the grid and then doing the integrals explicitly. 
The maximal value of M was 640. The achieved accuracy was about 10~ 5 . 

Thus calculated values of f n (y) were summed over n and integrated over v as indicated 
in (16)-(18) to obtain E and D. Stable results were obtained with the quite high maximal 
values n m = 300 and v m = 150.. In the part of n, v space outside the rectangle < n < 
30, < v < 15. the asymptotic expressions (22) and (23) were used for f n (v). 

As a result we calculated E and D as a quadratic form in the variational parameters 
Cfc. Afterwards the minimal value e of E subject to condition D = 1 was found by standard 
methods. 

The results for different number of parameters ./V are presented in the Table. To see the 
importance of the high n, v region, we also present the values of energy e± calculated without 
the asymptotical region, that is restricting to < n < 30, < v < 15.. The standard 
precision corresponds to the r,(f> grid 320x320. To clarify the accuracy achieved we also 
present the results with a double precision (the grid 640x640) for N = 5. 

Inspecting these results we see that the final accuracy in energy is of the order 7.10 -5 . This 
implies that taking N > 6 has no sense within the precision achieved, since the corresponding 
change in energy is of the same order or less. 
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So our conclusion is that the variational odderon energy with a trial function (13) is 
given by (5) and that with the accuracy achieved in the course of numerical integration, as 
described above, the maximal number of terms to be taken in the trial function is N = 6, 
although already with N = 3 used in [8] the energy is obtained up to 1%. 
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7 Table 

Odderon energy and parameters of the trial function 



N 


e 


C2,C 3 , ...(ci = 1) 


ei 


3 


0.22865 


-0.5036,0.2895 


0.21990 


4 


0.22632 


-0.2791,-0.3190,0.3609 


0.20663 


5 


0.22627 


-0.2009,-0.5052,0.1557, 
0.3779 


0.17919 


5* 


0.22634 


-0.2021,-0.5028,0.1543, 
0.3775 




6 


0.22619 


-0.3735,0.08842,-0.9765, 
1.003,0.3471 


-0.0015803 


7 


0.22618 




-0.37771 


8 


0.22616 




-0.54117 


9 


0.22616 




-0.60273 



*) double precision. 



